Evaluation of phylogenetic reconstruction methods using bacterial whole genomes: a simulation based study

Background: Phylogenetic reconstruction is a necessary first step in many analyses which use whole genome sequence data from bacterial populations. There are many available methods to infer phylogenies, and these have various advantages and disadvantages, but few unbiased comparisons of the range of approaches have been made. Methods: We simulated data from a defined "true tree" using a realistic evolutionary model. We built phylogenies from this data using a range of methods, and compared reconstructed trees to the true tree using two measures, noting the computational time needed for different phylogenetic reconstructions. We also used real data from Streptococcus pneumoniae alignments to compare individual core gene trees to a core genome tree. Results: We found that, as expected, maximum likelihood trees from good quality alignments were the most accurate, but also the most computationally intensive. Using less accurate phylogenetic reconstruction methods, we were able to obtain results of comparable accuracy; we found that approximate results can rapidly be obtained using genetic distance based methods. In real data we found that highly conserved core genes, such as those involved in translation, gave an inaccurate tree topology, whereas genes involved in recombination events gave inaccurate branch lengths. We also show a tree-of-trees, relating the results of different phylogenetic reconstructions to each other. Conclusions: We recommend three approaches, depending on requirements for accuracy and computational time. Quicker approaches that do not perform full maximum likelihood optimisation may be useful for many analyses requiring a phylogeny, as generating a high quality input alignment is likely to be the major limiting factor of accurate tree topology. We have publicly released our simulated data and code to enable further comparisons

Synergistic activity of mobile genetic element defences in Streptococcus pneumoniae

A diverse set of mobile genetic elements (MGEs) transmit between Streptococcus pneumoniae cells, but many isolates remain uninfected. The best-characterised defences against horizontal transmission of MGEs are restriction-modification systems (RMSs), of which there are two phase-variable examples in S. pneumoniae. Additionally, the transformation machinery has been proposed to limit vertical transmission of chromosomally integrated MGEs. This work describes how these mechanisms can act in concert. Experimental data demonstrate RMS phase variation occurs at a sub-maximal rate. Simulations suggest this may be optimal if MGEs are sometimes vertically inherited, as it reduces the probability that an infected cell will switch between RMS variants while the MGE is invading the population, and thereby undermine the restriction barrier. Such vertically inherited MGEs can be deleted by transformation. The lack of between-strain transformation hotspots at known prophage att sites suggests transformation cannot remove an MGE from a strain in which it is fixed. However, simulations confirmed that transformation was nevertheless effective at preventing the spread of MGEs into a previously uninfected cell population, if a recombination barrier existed between co-colonising strains. Further simulations combining these effects of phase variable RMSs and transformation found they synergistically inhibited MGEs spreading, through limiting both vertical and horizontal transmission

Variation at the capsule locus, cps, of mistyped and non-typable Streptococcus pneumoniae isolates

The capsule polysaccharide locus (cps) is the site of the capsule biosynthesis gene cluster in encapsulated Streptococcus pneumoniae. A set of pneumococcal samples and non-pneumococcal streptococci from Denmark, the Gambia, the Netherlands, Thailand, the UK and the USA were sequenced at the cps locus to elucidate serologically mistyped or non-typable isolates. We identified a novel serotype 33B/33C mosaic capsule cluster and previously unseen serotype 22F capsule genes, disrupted and deleted cps clusters, the presence of aliB and nspA genes that are unrelated to capsule production, and similar genes in the non-pneumococcal samples. These data provide greater understanding of diversity at a locus which is crucial to the antigenic diversity of the pathogen and current vaccine strategies

Dominant Role of Nucleotide Substitution in the Diversification of Serotype 3 Pneumococci over Decades and during a Single Infection

Streptococcus pneumoniae of serotype 3 possess a mucoid capsule and cause disease associated with high mortality rates relative to other pneumococci. Phylogenetic analysis of a complete reference genome and 81 draft sequences from clonal complex 180, the predominant serotype 3 clone in much of the world, found most sampled isolates belonged to a clade affected by few diversifying recombinations. However, other isolates indicate significant genetic variation has accumulated over the clonal complex’s entire history. Two closely related genomes, one from the blood and another from the cerebrospinal fluid, were obtained from a patient with meningitis. The pair differed in their behaviour in a mouse model of disease and in their susceptibility to antimicrobials, with at least some of these changes attributable to a mutation that upregulated the patAB efflux pump. This indicates clinically important phenotypic variation can accumulate rapidly through small alterations to the genotype

Meningococcal genetic variation mechanisms viewed through comparative analysis of Serogroup C strain FAM18

Copyright @ 2007 Public Library of ScienceThe bacterium Neisseria meningitidis is commonly found harmlessly colonising the mucosal surfaces of the human nasopharynx. Occasionally strains can invade host tissues causing septicaemia and meningitis, making the bacterium a major cause of morbidity and mortality in both the developed and developing world. The species is known to be diverse in many ways, as a product of its natural transformability and of a range of recombination and mutation-based systems. Previous work on pathogenic Neisseria has identified several mechanisms for the generation of diversity of surface structures, including phase variation based on slippage-like mechanisms and sequence conversion of expressed genes using information from silent loci. Comparison of the genome sequences of two N. meningitidis strains, serogroup B MC58 and serogroup A Z2491, suggested further mechanisms of variation, including C-terminal exchange in specific genes and enhanced localised recombination and variation related to repeat arrays. We have sequenced the genome of N. meningitidis strain FAM18, a representative of the ST-11/ET-37 complex, providing the first genome sequence for the disease-causing serogroup C meningococci; it has 1,976 predicted genes, of which 60 do not have orthologues in the previously sequenced serogroup A or B strains. Through genome comparison with Z2491 and MC58 we have further characterised specific mechanisms of genetic variation in N. meningitidis, describing specialised loci for generation of cell surface protein variants and measuring the association between noncoding repeat arrays and sequence variation in flanking genes. Here we provide a detailed view of novel genetic diversification mechanisms in N. meningitidis. Our analysis provides evidence for the hypothesis that the noncoding repeat arrays in neisserial genomes (neisserial intergenic mosaic elements) provide a crucial mechanism for the generation of surface antigen variants. Such variation will have an impact on the interaction with the host tissues, and understanding these mechanisms is important to aid our understanding of the intimate and complex relationship between the human nasopharynx and the meningococcus.This work was supported by the Wellcome Trust through the Beowulf Genomics Initiative

Pneumococcal Colonization and Virulence Factors Identified Via Experimental Evolution in Infection Models

Streptococcus pneumoniae is a commensal of the human nasopharynx and a major cause of respiratory and invasive disease. We examined adaptation and evolution of pneumococcus, within nasopharynx and lungs, in an experimental system where the selective pressures associated with transmission were removed. This was achieved by serial passage of pneumococci, separately, in mouse models of nasopharyngeal carriage or pneumonia. Passaged pneumococci became more effective colonizers of the respiratory tract and we observed several examples of potential parallel evolution. The cell wall-modifying glycosyltransferase LafA was under strong selection during lung passage, whereas the surface expressed pneumococcal vaccine antigen gene pvaA and the glycerol-3-phosphate dehydrogenase gene gpsA were frequent targets of mutation in nasopharynx-passaged pneumococci. These mutations were not identified in pneumococci that were separately evolved by serial passage on laboratory agar. We focused on gpsA, in which the same single nucleotide polymorphism arose in two independently evolved nasopharynx-passaged lineages. We describe a new role for this gene in nasopharyngeal carriage and show that the identified single nucleotide change confers resistance to oxidative stress and enhanced nasopharyngeal colonization potential. We demonstrate that polymorphisms in gpsA arise and are retained during human colonization. These findings highlight how within-host environmental conditions can determine trajectories of bacterial evolution. Relative invasiveness or attack rate of pneumococcal lineages may be defined by genes that make niche-specific contributions to bacterial fitness. Experimental evolution in animal infection models is a powerful tool to investigate the relative roles played by pathogen virulence and colonization factors within different host niches

Genomic analysis and comparison of two gonorrhoea outbreaks

© 2016 Didelot et al.Gonorrhea is a sexually transmitted disease causing growing concern, with a substantial increase in reported incidence over the past few years in the United Kingdom and rising levels of resistance to a wide range of antibiotics. Understanding its epidemiology is therefore of major biomedical importance, not only on a population scale but also at the level of direct transmission. However, the molecular typing techniques traditionally used for gonorrhea infections do not provide sufficient resolution to investigate such fine-scale patterns. Here we sequenced the genomes of 237 isolates from two local collections of isolates from Sheffield and London, each of which was resolved into a single type using traditional methods. The two data sets were selected to have different epidemiological properties: the Sheffield data were collected over 6 years from a predominantly heterosexual population, whereas the London data were gathered within half a year and strongly associated with men who have sex with men. Based on contact tracing information between individuals in Sheffield, we found that transmission is associated with a median time to most recent common ancestor of 3.4 months, with an upper bound of 8 months, which we used as a criterion to identify likely transmission links in both data sets. In London, we found that transmission happened predominantly between individuals of similar age, sexual orientation, and location and also with the same HIV serostatus, which may reflect serosorting and associated risk behaviors. Comparison of the two data sets suggests that the London epidemic involved about ten times more cases than the Sheffield outbreak. IMPORTANCE: The recent increases in gonorrhea incidence and antibiotic resistance are cause for public health concern. Successful intervention requires a better understanding of transmission patterns, which is not uncovered by traditional molecular epidemiology techniques. Here we studied two outbreaks that took place in Sheffield and London, United Kingdom. We show that whole-genome sequencing provides the resolution to investigate direct gonorrhea transmission between infected individuals. Combining genome sequencing with rich epidemiological information about infected individuals reveals the importance of several transmission routes and risk factors, which can be used to design better control measures

Microevolution of extensively drug-resistant tuberculosis in Russia.

Extensively drug-resistant (XDR) tuberculosis (TB), which is resistant to both first- and second-line antibiotics, is an escalating problem, particularly in the Russian Federation. Molecular fingerprinting of 2348 Mycobacterium tuberculosis isolates collected in Samara Oblast, Russia, revealed that 72%belonged to the Beijing lineage, a genotype associated with enhanced acquisition of drug resistance and increased virulence. Whole-genome sequencing of 34 Samaran isolates, plus 25 isolates representing global M. tuberculosis complex diversity, revealed that Beijing isolates originating in Eastern Europe formed a monophyletic group. Homoplasic polymorphisms within this clade were almost invariably associated with antibiotic resistance, indicating that the evolution of this population is primarily driven by drug therapy. Resistance genotypes showed a strong correlation with drug susceptibility phenotypes. A novel homoplasic mutation in rpoC, found only in isolates carrying a common rpoB rifampicin-resistance mutation, may play a role in fitness compensation. Most multidrug-resistant (MDR) isolates also had mutations in the promoter of a virulence gene, eis, which increase its expression and confer kanamycin resistance. Kanamycin therapy may thus select for mutants with increased virulence, helping preserve bacterial fitness and promoting transmission of drug-resistant TB strains. The East European clade was dominated by two MDR clusters, each disseminated across Samara. Polymorphisms conferring fluoroquinolone resistance were independently acquired multiple times within each cluster, indicating that XDR TB is currently not widely transmitted. © 2012 by Cold Spring Harbor Laboratory Press

Investigation of sequential outbreaks of Burkholderia cepacia and multidrug-resistant extended spectrum β-lactamase producing Klebsiella species in a West African tertiary hospital neonatal unit: a retrospective genomic analysis

Background Sick newborns admitted to neonatal units in low-resource settings are at an increased risk of developing hospital-acquired infections due to poor clinical care practices. Clusters of infection, due to the same species, with a consistent antibiotic resistance profile, and in the same ward over a short period of time might be indicative of an outbreak. We used whole-genome sequencing (WGS) to define the transmission pathways and characterise two distinct outbreaks of neonatal bacteraemia in a west African neonatal unit. Methods We studied two outbreaks of Burkholderia cepacia and multidrug-resistant extended spectrum β-lactamase (ESBL)-producing Klebsiella pneumoniae in a neonatal unit that provides non-intensive care on the neonatal ward in the Edward Francis Small Teaching Hospital, Banjul, The Gambia. We used WGS to validate and expand findings from the outbreak investigation. We retrospectively sequenced all clinical isolates associated with each outbreak, including isolates obtained from swabs of ward surfaces, environmental fluid cultures, intravenous fluids, and antibiotics administered to newborns. We also sequenced historical B cepacia isolates associated with neonatal sepsis in the same ward. Results Between March 1 and Dec 31, 2016, 321 blood cultures were done, of which 178 (55%) were positive with a clinically significant isolate. 49 episodes of neonatal B cepacia bacteraemia and 45 episodes of bacteraemia due to ESBL-producing K pneumoniae were reported. WGS revealed the suspected K pneumoniae outbreak to be contemporaneous outbreaks of K pneumoniae (ST39) and previously unreported Klebsiella quasipneumoniae subspecies similipneumoniae (ST1535). Genomic analysis showed near-identical strain clusters for each of the three outbreak pathogens, consistent with transmission within the neonatal ward from extrinsically contaminated in-use intravenous fluids and antibiotics. Time-dated phylogeny, including retrospective analysis of archived bacterial strains, suggest B cepacia has been endemic in the neonatal ward over several years, with the Klebsiella species a more recent introduction. Interpretation Our study highlights the emerging threat of previously unreported strains of multidrug-resistant Klebsiella species in this neonatal unit. Genome-based surveillance studies can improve identification of circulating pathogen strains, characterisation of antimicrobial resistance, and help understand probable infection acquisition routes during outbreaks in newborn units in low-resource settings. Our data provide evidence for the need to regularly monitor endemic transmission of bacteria within the hospital setting, identify the introduction of resistant strains from the community, and improve clinical practices to reduce or prevent the spread of infection and resistance